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ABSTRACT 

A summary of some recent results from lattice simulations is presented. These 
include first calculations of the strangeness magnetic moment of the nucleon, 
three new studies of the gluon propagator, flux tube attraction in J7(l) gauge 
theory, and the static-quark potential and its gluonic excitations. 



1. Introduction 

Lattice simulations afford a means of studying confinement from first principles. 
Faster computers and more efficient simulation techniques (duality transformations, 
improved actions, anisotropic lattices) are allowing access to physics in previously un- 
explored territories: the strangeness content of the nucleon, the infrared region of the 
Landau-gauge gluon propagator, the long flux in U(l), the rich glueball spectrum, and 
the static-quark potential and its gluon excitations for quark- ant iquark separations 
as large as 4 fm. In this talk, the new lattice results presented at this conference are 
summarized. 

2. Strangeness magnetic moment of the nucleon 

The strangeness content of the nucleon is currently of much interest. A significant 
ss content in the nucleon can resolve the long-standing discrepancy between the pion- 
nucleon sigma term extracted from low-energy pion-nucleon scattering and that from 
the octet baryon masses. A large and negative strange-quark polarization has emerged 
from spin structure function studies at CERN and SLAC, combined with neutron and 
hyperon (3 decays. Substantial contributions from strange quarks in the axial-vector 
and scalar current matrix elements would seem to imply similar such contributions to 
the vector current, but this would jeopardize the SU(6) prediction of the neutron to 
proton magnetic moment ratio, which lends credence to the valence quark picture. The 
SAMPLE collaboration has recently measured the neutral weak magnetic form factor 
in elastic parity-violating electron scattering at backward angles, and obtained the 
strangeness magnetic form factor by subtracting out the nucleon magnetic form factor. 
Elastic ep and e 4 He parity- violation experiments are currently planned at TJNAF to 
measure the asymmetry at forward angles to extract the strangeness electric mean- 
square radius. 

Williams presented resultstl from the first lattice calculation of the strangeness 
magnetic G s M (q 2 ) and electric G s E (q 2 ) form factors, where q is the four- momentum 
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Fig. 1. Magnetic G^ f (q 2 ) and electric G^(q 2 ) form factors for the proton (N = p) and the 
neutron (N = n), where q is the four- momentum transfer. The open circles indicate results 
from the connected insertion only, while the full results are shown as solid circles. The open 
circles are shifted slighted to the left to avoid overlap. The solid line is the fit to experiment. 



transfer. These form factors were obtained from ratios of two- and three-point cor- 
relation functions involving the vector current. The simulations were carried out on 
a 16 3 x 24 lattice with an inverse spacing a -1 = 1.72(4) GeV and were possible only 
because of recent developments in Z 2 noise techniques for quark propagators. The calcu- 
lations were done in the quenched approximation; contributions due to sea quarks were 
explicitly included at one fermion-loop order (the disconnected insertion diagram). The 
connected insertion diagram yields the contributions from the valence quarks and their 
Z graphs. Three values for the Wilson hopping parameter were used, corresponding to 
quark masses of 120, 200, and 360 MeV. 

Results for G s M (q 2 ) and G s E (q 2 ) were obtained at four nonzero values of q 2 . Ex- 



trapolating to the limit q 2 — >0 using a single pole form yielded G S M (0) = —0.36 ± 0.20 
and a value for G S E (0) consistent with zero. The electric mean-square radius was 
{ r s)E = —0.061 ± 0.003 fm 2 . For the u and d quarks in the disconnected insertion dia- 
gram, di s (0) = —0.65 ± 0.30 was found. Adding all of the sea quark contributions 
from the disconnected insertion diagram gives fidis = (2/3G^ dis (0) — l/3G^ dis (0) — 
1/3G s m (0))/jln = — 0.097±0.037/ijv, which, when combined with the connected insertion 
diagram, brings the ratio of neutron to proton magnetic moments \i n j fi p = — 0.68±0.04 
into agreement with the experimental value of 0.685. The accidental cancellation be- 
tween the disconnected insertion diagram and the rest-frame Z-graph component of 
the connected insertion diagram leaves the nonrelativistic valence contribution domi- 
nant, which explains why the SU (6) valence quark picture works well for the \i n f fi p 
ratio; SU(6) fails in the axial-vector and scalar cases because this cancellation does 
not occur. The results for the form factors are shown in Fig. |I[ Contributions from 
the strange quarks are not particularly small, but they are largely cancelled by the u 
and d quark contributions, resulting in a small overall effect to the form factors (the 
effect is greatest in the nucleon electric form factor). Issues related to the quenched 
approximation, chiral extrapolations, finite lattice-spacing artifacts, and finite volume 
errors must still be explored. 

3. Gluon propagator 

The connection between confinement and the gauge- dependent gluon propagator 
is unclear. However, knowledge of the infrared behaviour of the gluon propagator is 
important for various approaches to modelling confinement. Some studies based on 
Dyson-Schwinger equations claim that an infrared enhanced (such as 1/g 4 ) propagator 
is required for confinement; others assume that a dynamically generated gluon mass 
leads to an infrared finite behaviour. Lattice simulations can tell us the true infrared 
behaviour of the gluon propagator from first principles. However, large finite-volume 
effects have prevented previous lattice studies from accessing the momentum region 
where the relevant nonperturbative behaviour is expected. At this conference, three 
groups reported on new lattice studies of the gluon propagator; two worked in the 
Landau gauge and one used the maximally-abelian gauge in SU{2) with residual U(l) 
Landau gauge fixing. 

Williams presented results,^ shown in Fig. 0(a), revealing new structure in the 
gluon propagator for momenta as small as 0.4 GeV. Three separate simulations were 
done in order to control systematic errors from finite lattice spacing and finite volume. 
They obtained results in the quenched approximation on a large (32 3 x 64) lattice 
for a lattice spacing near 0.1 fm, which enabled them to probe momenta as small as 
400 MeV. Gribov copies were ignored. They also verified the tensor structure of the 
propagator, given in the continuum by 

DfM = 5 ah (<W " qj f) D(q% (1) 

where q is the four-momentum of the gluon, a and b are colour indices, and /i and v 
are the space-time indices. One of their main findings was to rule out a 1/g 4 infrared 
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Fig. 2. (a) The Landau-gauge gluon propagator in momentum space multiplied by q 2 , with 
nearby points averaged J3 The line is a best fit using the form given in Eq. (^). The scale is set 
using the string tension, (b) The SU{2) gluon propagator^ in the maximally-abelian gauge 
with residual U{\) Landau gauge fixing. The colour diagonal and off-diagonal components 
are shown. 



enhancement (for the case when quark-antiquark pair creation is neglected). They fit 
several phenomenological functions to their data, but most resulted in very large x 2 
values. In particular, the Gribov, Stingl, Marenzoni, and Cornwall forms were all shown 
to fail. One of the few functions to provide a satisfactory description of the data over 
a wide range of momenta was found empirically to be 



2 , nsr2 \ - d o s 



A 1/1, q 2 + M{j Y 



D(q) ~ Z + (Mi R ) 1+a + ^Mf w U ln ^-J I' (2) 

where Z = 1.781& A = 0.49^ 7 , aM uv = 0.20t?£, aM m = 0A3tl and a = 0.95tl 
were the fit parameters, and a -1 = 1.885 GeV and dp = 13/44. This function does 
incorporate some ultraviolet information from perturbation theory. In the future, these 
authors intend to use an improved-gauge field action and increase the lattice volume 
in order to probe deeper into the infrared. 

Results for the SU (2) gluon propagator in the maximally-abelian gauge (with 
residual U(l) Landau gauge fixing) were presented by SuganumaQ and are shown in 
Fig. ^|(b) . Two simulations on a 12 3 x 24 lattice were done, but the lattice spacings were 
too similar to allow any serious study of systematic errors from finite spacing and finite 
volume. Their focus was to demonstrate that the colour off-diagonal component has 
a Yukawa fall-off corresponding to a large mass so that long-range physics (r > 0.35 
fm) must be dominated by the colour diagonal component. This, they claimed, is 
the origin of abelian dominance in nonperturbative QCD. The dependence of the off- 
diagonal gluon's effective mass on the residual U(l) gauge fixing was not studied. It 
should be noted here that abelian dominance of nonperturbative QCD has not yet been 
convincingly established, and that confinement is much more than the area law of the 
Wilson loop. 
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Fig. 3. Results from dual simulations of flux tubes in U(l). (a) Ratio of the string tension 
of a doubly-charged toroidal flux tube over that of a singly-charged flux tube with respect 
to the coupling (5. (b) Longitudinal electric field profile of two interacting flux tubes in the 
symmetry plane (E\\, solid line) for (3 = 0.96. The equal charges are separated by 4a, where 
a is the lattice spacing. Dotted lines are field profiles E\\i and E^ 2 for single flux tubes at 
x = —2a and x = +2a. The noninteracting superposition E\\\ + Ep is shown as a dashed 
line. 



Nakajima and FuruiP presented preliminary results from not only quenched sim- 
ulations, but also from Langevin simulations including dynamical fermions (for three 
values of the Wilson hopping parameter k = 0.1, 0.15, and 0.2). Smeared gauge fixing 
was used to resolve Gribov ambiguities. Their lattices were very small (4 3 x 8) and 
much more work is needed to control systematic errors before any conclusions can be 
drawn on the role of quark-antiquark pair creation in the gluon propagator. 

4. Dual lattice simulations of U(l) flux tubes 

The study of confinement in the strong-coupling phase of compact U(l) in four 
dimensions can be viewed as a stepping stone on the way to understanding confinement 
in QCD. Confinement in the U(l) theory is driven by magnetic monopoles whose 
currents squeeze the electric field into flux tubes. A dual transformation of the £7(1) 
path integral exists so that simulations can be carried out much more efficiently in 
the dual theory, which corresponds to a certain limit of a dual Higgs model. Results 
obtained in the dually-transformed theory were presented^ at this conference, along 
with some first computations in the dual Abelian Higgs model. This study focussed on 
the question of whether the vacuum corresponds to a type-I or a type-II superconductor, 
and the role of quantum fluctuations of the dual degrees of freedom were investigated. 

First, it was shown that the string tension scales proportionally to the charge 
rather than the square of the charge in the confinement phase. This is shown in Fig. |3](a) 
in which the string tension of a closed double flux tube is compared to that of a single 
closed flux tube as a function of the coupling (3. These results were obtained on a 
8 3 x 16 lattice. In the Coulomb phase (large (3 values), the ratio of the string tension 
of the double flux tube to that of the single flux tube is 4 in agreement with the 
expected result for the energy of a homogeneous field of double strength; below the 
transition in the confinement phase, the ratio is 2, in agreement with the expectations 



for a dual superconductor. Flux tube attraction in the vicinity of the phase transition 
was confirmed by an examination of the longitudinal electric field profile at f3 — 0.96 
as shown in Fig. 0(b). Two parallel flux tubes of length 22a, where a is the lattice 
spacing, were brought next to each other, a transverse distance 4a apart, and the 
resulting longitudinal electric field in the symmetry plane was measured. The profile 
differed dramatically from the noninteracting superposition of two single flux tubes; the 
formation of one flux tube was observed, in analogy to a type-I superconductor. These 
authors also compared their simulation results to the predictions from a dual London 
equation. Agreement was found for small charge distances, but not for distances greater 
than ten lattice spacings (they probed lengths up to 20a). They interpreted this as 
further evidence that the confining U(l) behaves as an effective type-I superconductor. 

5. Static-quark potentials 

Three groups reported on calculations of static quark- ant iquark potentials; four- 
quark energies in SU (2) lattice gauge theory were also presented. 

The static-quark potential and its gluon excitations are very useful probes of 
confinement. It is generally believed that at large quark-antiquark separation r, the 
linearly-growing ground-state energy of the glue is the manifestation of the confining 
flux whose fluctuations can be described in terms of an effective string theory. The 
lowest-lying excitations are then the Goldstone modes associated with spontaneously- 
broken transverse translational symmetry. Expectations are less clear for small r. The 
gluon excitations of the static-quark potential are also useful for studying hybrid heavy- 
quark mesons using a Born-Oppenheimer expansion. 

The first comprehensive determination of the low-lying spectrum of gluonic ex- 
citations in the presence of a static quark-antiquark pair was presented at this con- 
ference B The glue energies for r from 0.1 to 4 fm were extracted from Monte Carlo 
estimates of generalized Wilson loops in eight simulations using an improved gauge- 
field action. The use of anisotropic lattices in which the temporal lattice spacing at was 
much smaller than the spatial spacing a s was crucial for resolving the glue spectrum, 
particularly for large r. Finite volume errors were shown to be negligible. Particular at- 
tention was paid to the volumes for very large r: not only were checks carried out using 
additional simulations, but also by studying volume effects for a naive Nambu-Goto 
string. Agreement of energies obtained using different quark-antiquark orientations on 
the lattice was used to check the smallness of finite-spacing errors and to help identify 
the continuum rotational quantum numbers corresponding to each level (there are only 
discrete symmetries on the lattice). Results for lattice spacings ranging from 0.12 to 
0.29 fm were obtained to facilitate extrapolation to the continuum limit. The hadronic 
scale parameter r ~ 0.5 fm was used to set the scale. 

The continuum-limit extrapolations are shown in Fig. |j. The energies are labelled 
by the magnitude A of the projection of the total angular momentum of the gluons 
onto the molecular axis, and the symmetry t]cp under spatial inversion with charge 
conjugation. The capital Greek letters S, IT, A, . . . are used to indicate states with 
A = 0, 1, 2, 3, . . ., respectively. States with r\ C p = 1(— 1) are denoted by the subscripts g 
(u). The S states are additionally labelled by a superscript + (— ) depending on whether 
they are even (odd) under a reflection in a plane containing the molecular axis. The 
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Fig. 4. Gluon excitations of the static-quark potential. The dashed lines indicate the locations 
of the rrnr/r gaps above the curve for m = 1, 2, 3, and 4. The dotted curves are the naive 
Nambu-Goto energies in four-dimensions. Note that a Q u interpretation for the curve labelled 
H' u cannot be ruled out. 

ground-state E+ is the familiar static-quark potential. The lowest-lying excitation is 
the Ii u . There is definite evidence of a band structure at large r: the H' , U g , and A g 
form the first band above the U u ; the E+, £~, U' u /§ u , and A u form another band; and 
the S~ is the highest level. The level orderings and approximate degeneracies of the 
gluon energies at large r match, without exception, those expected of the Goldstone 
excitations which are a universal feature of any effective string theory description of 
the long confining flux. However, the precise rrnr/r gap behaviour is not observed. For 
separations less than 2 fm, the gluon energies lie well below the Goldstone energies and 
the Goldstone degeneracies are no longer observed. The two X - states are in violent 
disagreement with expectations from a fluctuating string. These somewhat surprising 
results cast serious doubts on the validity of treating glue in terms of a fluctuating 
string for quark-antiguark separations less than 2 fm. 

Some authorsQ'il claim that the glue can be described merely as a mathematical 
string using the Nambu-Goto action (turning a blind eye to the fact that the Nambu- 
Goto string cannot be consistently quantized in four dimensions). These results clearly 
contradict this claim. For r greater than 1 fm, the Il u and IT 9 energies disagree the 
least with the naive Nambu-Goto energies. In all other cases, the discrepancies are 
significant. The E~ and S~ levels look nothing like the predictions from the naive 
Nambu-Goto string. The departures of these levels from the Nambu-Goto energies are 
so severe that explanations in terms of mixings with other string modes or consequences 



of small-r symmetry requirements are difficult to accept. Rather, such behaviour signals 
the complete failure of the Nambu-Goto intrepretation. 

Bakers presented a comparison of results on the central, spin- and momentum- 
dependent interquark potentials from a dual superconductor model with those deter- 
mined from first principles using lattice simulations. The model is an effective theory 
of long-distance Yang-Mills in which an octet of dual potentials are coupled minimally 
to three octets of scalar Higgs fields carrying colour magnetic charge. Monopole con- 
densation takes place and the dual potentials couple to a quark-antiquark pair via a 
Dirac string connecting the pair. Chromoelectric flux is confined as the quark-antiquark 
separation increases as a result of the dual Meissner effect. Assuming a particular spon- 
taneous symmetry breaking sequence, the Lagrangian of the model can be replaced by 
an Abelian Higgs Lagrangian. The model has two parameters, a coupling constant a s 
and the vacuum expectation value of the dual Higgs field 0o- The interquark poten- 
tials were calculated from the Wilson loop of the dual theory. The two parameters of 
the model were determined by fitting to the lattice results for the central potential. 
The remaining nine potentials are then uniquely determined and were compared to the 
lattice data. They were shown to be in fair agreement for r greater than 0.2 fm; the 
agreement was aided by large uncertainties in the simulation results. 

KlindworthHI calculated the static quark potential in transverse light-front QCD. 
In this approach, the two light-cone coordinates x = (x° are continuous 

and their associated degrees of freedom are non-compact, while the other two trans- 
verse directions x± are discretized and the gauge fields associated with these directions 
are expressed in terms of compact, non-unitary link variables. The calculation was 
carried out in the large N c limit, where iV c is the number of colours. In so doing, 
quark-antiquark pair creation could be ignored. The static quark potential was then 
calculated using Lanczos diagonalization techniques which necessitated a Fock-space 
truncation. A linearly confining potential was found, and an excited potential was also 
determined. The calculation served mainly as an initial test of this new regularization 
and computational scheme. 

Furuillll presented results for four-quark energies from quenched SU(2) simula- 
tions. The aim of the analysis was to obtain a compact expression for the four-quark 
interactions for an arbitrary spatial configuration of the quarks. They have examined 
square, rectangular, tilted rectangular, linear, quadrilateral, non-planar (at least one 
link on axis), and tetrahedral (two links off axis) geometries on a 16 3 x 32 lattice at 
two different lattice spacings. Their results suffer from significant lattice artifacts. 

6. Conclusion 

A plethora of new results from lattice simulations are shedding light in some 
of the previously dark corridors of the labyrinth surrounding confinement. Computer 
simulations of gluons and quarks remain an important tool for navigating this maze. 
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